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Abstract 



We analyze the indirect detection of neutralino dark matter in supergravity scenar- 
ios with non-universal soft scalar and gaugino masses. In particular, the gamma-ray 
Q_i! flux arising from the galactic center due to neutralino annihilation is computed. In 

important regions of the parameter space it can be increased significantly (about two 
orders of magnitude) with respect to the universal scenario. This result is compared 
with the sensitivity of current and planned experiments, such as the satellite-based 
detectors EGRET and GLAST and the ground-based atmospheric telescope HESS. 
For example, for tan (3 < 35 the sensitivity region of GLAST and HESS can only be 
reached in non-universal scenarios. 
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1 Introduction 



There are substantial evidences suggesting the existence of dark matter in the Uni- 
verse. In addition to the observations of galactic rotation curves [1], cluster of galaxies, 
and large scale flows [2], implying 0.1 < QuMh 2 < 0.3, the recent data obtained by the 
WMAP satellite [3] confirm that dark matter must be present with 0.094 < Q,DMh 2 < 0.129. 

The lightest neutralino, Xi> a particle predicted by the supersymmetric (SUSY) 
extension of the standard model, is one of the leading candidates for dark matter. It is 
specially interesting because, in addition to be usually stable, it is a weakly interacting 
massive particle (WIMP), and therefore it can be present in the right amount to explain 
the observed matter density. 

There are promising methods for the indirect detection of WIMPs through the anal- 
ysis of their annihilation products [4] . One of them consists of detecting the gamma rays 
produced by these annihilations in the galactic halo [5]- [24]. For this one uses atmo- 
spheric Cherenkov telescopes or space-based gamma-ray detectors. In fact, one of the 
space-based experiments, the Energetic Gamma-Ray Experiment Telescope (EGRET) 
on the Compton Gamma-Ray Observatory detected a signal [25] that apparently can- 
not be explained with the usual gamma-ray background. In particular, after five years 
of mapping the gamma-ray sky up to an energy of about 20 GeV, there are evidences 
of a source of gamma rays above about 1 GeV in the galactic center region with a value 
for the flux of about 10 -8 cm' 2 s -1 . It is worth noticing, however, that alternative 
explanations for this result have been proposed modifying conveniently the standard 
theory of galactic gamma ray [26] . 

Projected experiments might clarify the situation. For example, a larger sensitivity 
will be provided by the Gamma-ray Large Area Space Telescope (GLAST) [27], which 
is scheduled for launch in 2006. In addition, energies as high as 300 GeV can be 
measured testing models with heavy dark matter more easily than EGRET. Even 
if alternative explanations solve the data problem pointed out by EGRET, GLAST 
will be able to detect a smaller flux of gamma rays from dark matter ~ 10~ 10 cm -2 
s _1 . Besides, ground-based atmospheric Cherenkov telescopes are complementary to 
space-based detectors. In particular, the High Energy Stereoscopic System (HESS) 
experiment [28], which has begun operations with four telescopes, has a large effective 
area of 10 5 m 2 (versus 1 m 2 for GLAST) allowing the detection of very high energy 
fluxes > 1 TeV. On the contrary, the energy threshold is larger, > 60 GeV. 

Given the experimental situation described above, concerning the detection of dark 
matter, and assuming that this is made of neutralinos, it is natural to wonder how big 
the annihilation cross section for its indirect detection can be. Obviously, this analysis 
is crucial in order to know the possibility of detecting dark matter in the experiments. 

This has been carried out [4, 16, 17, 18, 19] in the usual minimal supergravity 
(mSUGRA) scenario, where the soft terms of the minimal supersymmetric standard 
model (MSSM) are assumed to be universal at the unification scale, Mqut ~ 2 x 10 16 
GeV, and radiative electroweak symmetry breaking is imposed. Neutralino annihila- 
tions in the galactic center region seem not to be able to explain the excess detected 
by EGRET. In fact, only small regions of the parameter space of mSUGRA may be 
compatible with the sensitivity of GLAST and HESS, for standard dark matter den- 
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sity profiles. Clearly, in this case, more sensitive detectors producing further data 
are needed. However, this result might be modified by taking into account possible 
departures from the mSUGRA scenario. In particular, a more general situation in 
the context of SUGRA than universality, the presence of non-universal soft scalar and 
gaugino masses, might increase the annihilation cross section significantly, producing 
larger gamma-ray fluxes. 

The aim of this paper is to investigate this general case, where non-universalities 
are present in the scalar and gaugino sectors 1 , and to carry out a detailed analysis of 
the prospects for the indirect detection of neutralino dark matter in these scenarios. In 
the light of the recent experimental results, we will be specially interested in studying 
how big the gamma-ray flux from the galactic center can be. Our purpose is to provide 
a general analysis which can be used in the study of any concrete model. 

In this sense, it is worth noticing that such a non-universal structure can be re- 
covered in the low-energy limit of some phenomenologically appealing string scenarios. 
This is the case, for example, of heterotic orbifold models [29] or D-brane constructions 
in the type I string [30] . An analysis of the gamma-ray fluxes in these scenarios might 
be very interesting [31]. On the other hand, AMSB scenarios in an effective heterotic 
framework were studied in Ref. [24, 23]. 

The paper is organized as follows. In Section 2 we will discuss in general the 
gamma-ray detection from dark matter annihilation, paying special attention to its 
halo model dependence. The particle physics dependence will be reviewed in Section 3 
for the case of neutralinos in mSUGRA. In Section 4 we will study the general case 
where scalar and gaugino non-universalities are present. In particular, we will indicate 
the conditions under which a significant enhancement of the resulting gamma-ray flux 
is obtained. Finally, the conclusions are left for Section 5. 

2 Gamma-ray flux from dark matter annihila- 
tion 

Annihilation of dark matter particles in the galactic center [6], in the galactic halo [7] 
or in the haloes of nearby galaxies [8, 9] could produce detectable fluxes of high-energy 
photons. We are interested in the analysis of fluxes from the galactic center because 
they can be significantly enhanced for some specific dark matter density profiles [10], 
as we will see below. 

There are two possible types of gamma rays that can be produced by the annihi- 
lation. First, gamma-ray lines from processes XX ~^ 77 [H] an d XX ~^ 1% [12]. This 
signal would be very clear since the photons are basically mono-energetic. Unfortu- 
nately, the neutralino does not couple directly to the photon, the Feynman diagrams are 
loop suppressed, and therefore the flux would be small [13]. For more recent analyses of 
this possibility see Refs. [10, 18]. On the other hand, continuum gamma rays produced 
by the decay of neutral pions generated in the cascading of annihilation products will 
give rise to larger fluxes [14]. Although we will concentrate on the latter possibility in 

1 For analyses in the context of the eSMSSM scenario, where the parameters are defined directly at the 
electroweak scale, see e.g. the recent works [20, 21]. 
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the following, let us remark that we have also checked for all cases studied in the paper 
the flux of monochromatic gamma rays. This turns out to be too small for a standard 
NFW density profile, with an upper bound of about 10 -13 cm~ 2 s _1 , and therefore 
much below the sensitivity of GLAST. 

For the continuum of gamma rays, the differential flux coming from a direction 
forming an angle ip with respect to the galactic center is 

d$ 7 ^ 1 dN, 



dUdE 



= y l^L {(Tiv) _}_ f p 2 d i (2 .i) 



where the discrete sum is over all dark matter annihilation channels, dN^/dE^ is the 
differential gamma-ray yield, (<7ji>) is the annihilation cross section averaged over its 
velocity distribution, m x is the mass of the dark matter particle, and p is the dark 
matter density. It is worth noticing here that we have included in the above equation 
the factor 1/2 omitted in previous literature 2 . Assuming a spherical halo, p = p(r) with 
the galactocentric distance r 2 = I 2 + Rq — 21R$ cos ■0, where Rq is the solar distance to 
the galactic center (~ 8 kpc [32]). Following Ref. [10], one can separate in the above 
equation the particle physics part from the halo model dependence introducing the 
(dimensionless) quantity 

1 / 1 



J ^ = ^r( n.rv/ 3 / P \r(l^))dl. (2.2) 

8.5 kpc V0.3 GeV/cnrV J line of sight 

Thus the gamma-ray flux can be expressed as 

® 7 (E thr ) = 1-87 x 1(T 13 cm" 2 s" 1 

i E r^p. ( ("»>, \ fi«>«vy J(An)An (2 . 3 , 

o^L 7 d£L V lO-^cmSs" 1 / V m v / v ' y ' 



x 

' W lE thr 



Here Ethr is the lower threshold energy of the detector, and with respect to the upper 
limit of the integral notice that neutralinos move at galactic velocity and therefore 
their annihilation occurs at rest. The quantity J (AST) is defined as 



'(AO) = [ M ■ (2-4) 



Although the flux is maximized in the direction of the galactic center (corresponding 
to 4> = 0), rather than J(0) one must consider the integral of J(Y') over the spherical 
region of solid angle Afl given by the angular acceptance of the detector which is 
pointing towards the galactic center. Typically AQ is about 10~ 3 sr. 



2 As discussed in Ref. [9] the simplest way to see this is realizing that a is the annihilation cross section 
for a given pair of particles. Thus in a given volume with N WIMPs the annihilation rate will be av times 
the number of pairs N(N — l)/2. In the continuum limit this result can be approximated as avn 2 /2. 
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a (kpc) J(1(T 2 ) J(1(T 3 ) J(1(T 5 ) 


Isothermal 
Kravtsov et al. 

NFW 
Navarro et al. 
Moore et al. 


3.5 2.439 xlO 1 2.466 xlO 1 2.470 xlO 1 
10 1.625 xlO 1 1.932 xlO 1 2.368 xlO 1 
20 3.665 xlO 2 1.223 xlO 3 1.261 xlO 4 
0.952 xlO 3 3.015 xlO 3 1.407 xlO 4 
28 2.069 xlO 4 1.443 xlO 5 1.033 x 10 7 



Table 1: Different dark matter density profiles with the corresponding value of J (AVI), for 

An = icr 2 , icr 3 , icr 5 sr. 



2.1 Astrophysics contribution 

A crucial ingredient for the calculation of J, and therefore of the flux of gamma rays, 
is the dark matter density profile of our galaxy. The different profiles that have been 
proposed in the literature can be parameterised as [33] 

P{ ) (r/Rop[l + (r/a)°](/»-T)/° ' 1 ' 

where po is the local (solar neighborhood) halo density and a is a characteristic length. 
Although we will use po = 0.3 GeV/cm 3 throughout the paper, since this is just a 
scaling factor in the analysis, modifications to its value can be straightforwardly taken 
into account in the results. Obviously, for r = Rq one recovers the local density in 
Eq. (2.5). 

For (a, (3, j)=(2, 2, 0) one obtains the simple Isothermal halo model p{r) = po a 2 +r z ■, 
where a is the core radius. This profile is clearly non-singular at the galactic center 
r = 0. However, highly cusped profiles seems to be deduced from N-body simulations 
3 . Navarro, Frenk and White (NFW) [35] obtained (a, j3, 7)=(1, 3, 1), producing a 
profile with a behaviour 1/r at small distances. A more singular behaviour, 1/r 1 ' 5 , 
was obtained by Moore et al. [36] with (a, j3, 7)=(1.5, 3, 1.5). On the other hand, 
Kravtsov et al. [37] obtained a mild singularity towards the galactic center since for 
them (q,/3,7)=(2, 3,0.4). Moreover, Navarro et al. [38] have recently proposed a 
non-singular distribution 



P(r) = P-2 exp <! -- 



Q: 



r 

r-2 



1 



(2.6) 



with typical values a = 0.142, r_ 2 = 26.4 kpc, p_ 2 = 0.035 GeV cnr 3 . 

Clearly, the situation concerning the halo models is still unclear. The situation 
is similar from the observational viewpoint. Reconstructions of the observed rotation 
curves are claimed to be consistent or inconsistent with the predictions of N-body 
simulations depending on the authors. 

Taking all the above into account, we show in Table 1 the value of J(A$7), with 
AQ = 10~ 2 , 10~ 3 , 10 -5 sr, for the most widely used dark matter profiles. In the 
computation one usually assumes a constant density for r < 0.01 pc, in order to avoid 



3 For analytical derivations see e.g. the recent work [34], and references therein. 
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Figure 1: Differential photon yield dN 7 /dE 7 for several annihilation channels, using a WIMP with mass 
m w = 200 GeV. The total photon yield for each channel is also shown in parenthesis. 



the problem of the divergence of the singular profiles. Obviously, the gamma-ray flux 
can be largely enhanced depending on the profile chosen 4 . For example, for AO = 10 -3 
sr, the flux produced by the Moore et al. profile is about two orders of magnitude larger 
than the one produced by the NFW, and this is about two orders of magnitude larger 
than the flux produced by the Isothermal and Kravtsov profiles. 

For the sake of definiteness, we will use in our analysis a NFW profile with AO = 
10~ 3 sr. Using Table 1, the results can easily be extended to other profiles and angular 
resolutions, since J(AO) enters just as a scaling factor in the gamma-ray flux (2.3). 

2.2 Particle physics contribution 

Concerning the particle physics contribution to the gamma-ray flux (2.3), one has all 
the other factors in front of J(AO). Their values will depend on the specific theoretical 
scenario where the WIMPs are embedded. With respect to the differential photon 
yield, it is known that this shows scaling properties [14], i.e. it depends only on 
the scaled variable x = E^jm w , where m w is the mass of the WIMP. Thus can be 
reasonably parameterised for each possible annihilation channel of the WIMPs [14, 10]. 
For example, assuming that WIMPs annihilate producing WW, ZZ, bb, it, and uu, 
one can write dN^/dx = ae~ bx /x 1,5 , with (a, b) = (0.73,7.76) for the annihilation 

4 We are conservative and do not consider here the effect on this study of the possible existence of a super 
massive black hole at the center of our galaxy of 3.7 x 10 6 solar masses [39]. This might produce a 'spike' in 
the dark matter density profile leading to behaviour 1/r 2 4 [40], and therefore increasing substantially the 
gamma-ray flux [41]. 
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channels WW and ZZ [10, 42], (1.0,10.7) for bb, (1.1,15.1) for tt, and (0.95,6.5) for 
uu [4]. In fact, the spectrum of gamma-rays is very close for all the above channels 
as illustrated in Fig. 1, except in the case t + t~ which produces smaller values (see 
e.g. [14, 16]). In our computation of the gamma-ray flux below we use the more exact 
PYTHIA fragmentation Monte Carlo [43]. 

On the other hand, the influence of the theoretical framework is crucial for the value 
of the annihilation cross section o"j. As mentioned in the Introduction, we will work 
in the context of the simplest SUSY construction, where the lightest supersymmetric 
particle (LSP) is absolutely stable, and therefore a candidate for dark matter. It is 
remarkable that in most of the parameter space of this construction the LSP is an 
electrically neutral (also with no strong interactions) particle, called neutralino. This 
is welcome since otherwise the LSP would bind to nuclei and would be excluded as a 
candidate for dark matter from unsuccessful searches for exotic heavy isotopes. 

As a matter of fact, there are four neutralinos, x® (i = 1,2,3,4), since they are 
the physical superpositions of the fermionic partners of the neutral electroweak gauge 
bosons, called bino (B ) and wino (W3), and of the fermionic partners of the neutral 
Higgs bosons, called Higgsinos (H®, H%). Therefore the lightest neutralino, Xi, wm 
be the dark matter candidate. In the basis (B°, W®, H®, H%), the neutralino mass 
matrix is given by 

/ Mi — mz cos P sin 9 W mz sin (3 sin 9 W \ 

M2 mz cos (3 cos W —mz sin (3 cos 9 W 

— mz cos (3 sin 9 W mz cos (3 cos 9 W —fi 

\ mz sin (3 sin 9 W — mz sin (3 cos 9 W —fj, / 

(2-7) 

where Mi and M2 are the bino and wino masses respectively, ji is the Higgsino mass 
parameter and tan/3 = (H®)/(H^) is the ratio of Higgs vacuum expectation values. 
This can be diagonalized by a matrix Z such that we can express the lightest neutralino 
as 

X\ = Z U B° + Z 12 W° + Z 13 H° + Z lA H° . (2.8) 

It is commonly defined that Xi is mostly gaugino-like if P = \Zu\ 2 + \Z\2\ 2 > 0.9, 
Higgsino-like if P < 0.1, and mixed otherwise. 

In Fig. 2 we show the relevant Feynman diagrams contributing to neutralino annihi- 
lation. As we will see, the cross section can be significantly enhanced depending on the 
SUSY model under consideration. We will concentrate here on the SUGRA scenario, 
where the soft terms are determined at the unification scale, Mqut ~ 2 x 10 16 GeV, 
after SUSY breaking, and radiative electroweak symmetry breaking is imposed. 

Let us start reviewing the situation in the case of mSUGRA, where the soft terms 
are assumed to be universal. 

3 mSUGRA predictions for the gamma-ray flux 

Let us first recall that in mSUGRA one has only four free parameters defined at the 
GUT scale: the soft scalar mass m, the soft gaugino mass M, the soft trilinear coupling 
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Figure 2: Dominant ncutralino annihilation diagrams. Relevant parts of the amplitudes are shown explicitly 
Terms between parenthesis correspond to /„ and Z final states in second and fourth diagrams. V and N are 
chargino and neutralino mixing matrices. 



A, and tan 3. In addition, the sign of the Higgsino mass parameter, fi, remains also 
undetermined by the minimization of the Higgs potential, which at tree level implies 

„ m?r, — m 2 TT tan 2 3 1 

As is well known, in this scenario the lightest neutralino is mainly bino [44, 45]. 
To understand this result qualitatively recall the evolution of mr H towards generically 
large and negative values with the scale. Since /U 2 given by Eq. (3.9), for reasonable 
values of tan/?, can be approximated as, 

li 2 « -m 2 Hu - \m 2 , (3.10) 

then it becomes also large. In particular, becomes much larger than M\ and M2. 
Thus, as can be easily understood from Eqs. (2.7) and (2.8), the lightest neutralino will 
be mainly gaugino, and in particular bino, since at low energy M\ = | tan 2 6w^2 ~ 
0.5M 2 . 

Now, using Eq. (2.3) one can compute the gamma-ray flux for different values of the 
parameters. As a consequence of the Xi being mainly bino, the predicted Oi is small, 
and therefore the flux is below the present accessible experimental regions. Note that 
only Z\\ is large and then the contribution of diagrams in Fig. 2 will be generically 
small. In addition, the (tree-level) mass of the CP-odd Higgs, A, 

m\ = m 2 Hd +m 2 Hu +2n 2 , (3.11) 

will be large. The previous formula can be rewritten as 

™ 2 A ~m 2 Hd -m 2 Hu -M 2 , (3.12) 

where Eq. (3.10) has been used. As mentioned above, m 2 H evolves at low energy 
towards large and negative values, producing a large value for m^, and therefore a 
further suppression in the annihilation channels. 

This fact is shown in Fig. 3 using the DarkSUSY program [46] 5 . There, contours 
of the gamma-ray flux $ 7 in the parameter space (m, M) for tan/3 = 5, A = and 



5 As discussed in Eq. (2.1), we arc including in the computation the factor 1/2 which was missing in 
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H > are plotted, using the complete one-loop corrected effective potential. We will not 
consider in the calculation the opposite sign of ji because this would produce a negative 
contribution for the — 2, and, as will be discussed below, we are mainly interested 
in positive contributions. Recall that the sign of the contribution is basically given by 
fiM2, and that M, and therefore M2, can always be made positive after performing an 
U(1)r rotation. As input for the top mass, rntipole) = 175 GeV has been used. 

In the computation we are using a NFW profile with AQ = 1CT 3 sr, and recent 
experimental and astrophysical constraints are taken into account. In particular, the 
lower bounds on the masses of the supersymmetric particles and on the lightest Higgs 
have been implemented, as well as the experimental bounds on the branching ratio 
of the b — > 57 process and on a^ USY . Due to its relevance, the effect of the WMAP 
constraint on the dark matter relic density is shown explicitly. 

Concerning a^ USY , we have taken into account the recent experimental result for 
the muon anomalous magnetic moment [50], as well as the most recent theoretical 
evaluations of the Standard Model contributions [51]. It is found that when e + e~ 
data are used the experimental excess in (g M — 2) would constrain a possible super- 
symmetric contribution to be a^ USY = (27.1 ± 10) x 10~ 10 . At 2a level this implies 
7.1 x 10~ 10 < a^ USY < 47.1 x 10~ 10 . It is worth noticing here that when tau data are 
used a smaller discrepancy with the experimental measurement is found. In order not 
to exclude the latter possibility we will discuss the relevant value a^ USY = 7.1 x 10~ 10 
in the figures throughout the paper. The region of the parameter space to the right of 
this line will be forbidden (allowed) if one consider electron (tau) data. 

On the other hand, the measurements of B — > X s ~j decays at CLEO [52] and BELLE 
[53], lead to bounds on the branching ratio b — > 57. In particular we impose on our 
computation 2.33 x 10~ 4 < BR(b — > sj) < 4.15 x 10~ 4 , where the evaluation is carried 
out using the routine provided by the program micrOMEGAs [48]. This program is 
also used for our evaluation of a^ USY and relic neutralino density. 

As we can see in Fig. 3, the points corresponding to a sufficiently low relic neutralino 
density, and satisfying the Higgs mass bound 6 in mSUGRA, > 114 GeV, are located 
in the narrow coannihilation branch of the parameter space, i.e. the region where the 
stau is the next to the LSP producing efficient coannihilations. For these points the 
neutralino mass is larger than the top mass (note that we can deduce its value in the 
plot from the value of M, since in mSUGRA m^o ~ M\ ~ 0.4 M implying m^o > 200 
GeV), and tt final states are allowed. In fact, through the A-exchange channel, these are 
the main source of the gamma-ray flux. The contribution from bb states is also relevant, 
and the relation between both contributions can be obtained taking into account that 



previous literature, and also in the equations of the program. We have used the Fortran code SuSpect2 [47] 
to solve the renormalization group equations (RGE) for the soft supersymmetry breaking terms between the 
GUT scale and the electroweak scale. These parameters are then passed on to the C code micrOMEGAs 
[48] to perform the calculation of physical masses for the superpartners and various indirect constraints to 
be described below. We have compared our results with the ones obtained with the recent released version 
of DarkSUSY [49] , and they are very similar. Only for large values of tan (3 and M there is a factor ~ 3 of 
discrepancy in the flux. 

6 Let us remark that this bound must be relaxed accordingly when departures from mSUGRA are 
considered. 



8 



(see Fig. 2) (JA(tt) oc f3 t m 2 / tan 2 /?, <ta(M>) oc /3f,m 2 tan 2 /3, with /?/ = -m 2 /m|, 
and therefore <TA(it)/cyi(66) ~ t^^ ( rn t/ m b) 2 f° r sufficiently large neutralino mass. 
In particular, for tan/3 = 5 one obtains a A{tt) / a A{bb) ~ 2. There is of course the 
possibility of rr final states, but their contributions are negligible compared with the bb 
ones. Note that, following the above arguments, a a{tt) / a A{bb) ~ |(m r /mb) 2 ~ 0.05, 
where 3 is a colour factor. 

Let us finally remark that if we impose consistency with the e + e~ data concerning 
(g^ — 2), the whole parameter space would be excluded (notice that the region to the 
right of the black dashed line in Fig. 3 corresponds to a^ USY < 7.1 x 10 -10 ). 

We summarize the above results in Fig. 4, using the same parameter space as in 
Fig. 3. There, the values of <I> 7 allowed by all experimental constraints as a function 
of the neutralino mass m^o are depicted. Points shown with black stars have 0.094 < 
tt^oh 2 < 0.129, and are therefore favoured by WMAP. Points with 0.129 < U^oh 2 < 0.3 
are shown with light grey (magenta) triangles. We also consider the possibility that 
not all the dark matter is made of neutralinos, allowing 0.03 < ^^oh 2 < 0.094. In this 
case we have to rescale the density of neutralinos in the galaxy p(r) in Eqs. (2.5) and 
(2.6) by a factor O^o/i 2 /0.094. Points corresponding to this possibility are shown with 
dark grey (blue) boxes. 

We observe that, generically, the gamma-ray flux for a threshold of 1 GeV (left 
frame) is constrained to be <I> 7 < 10~ n cm~ 2 s _1 . Obviously, in this mSUGRA scenario 
with tan (3 = 5 the EGRET data cannot be reproduced, and more sensitive detectors 
producing further data are needed. Note that GLAST experiment, where values of 
the gamma-ray flux as low as 10 -10 cm~ 2 s" 1 will be accessible, is not sufficient. For 
atmospheric telescopes, such as HESS, the situation is similar, as can be seen in the 
plot on the right frame. 

The annihilation cross section can be increased, and therefore the gamma-ray flux, 
when the value of tan (3 is increased. Notice for instance that the contribution of 
the bb final state, which is proportional to tan 2 /3 is now more important than the 
tt one. For example, for tan (3 = 35, using the above approximation one obtains 
a A(ti) I 'a A(bb) ~ ^(rat/mb) 2 ~ 10~ 3 . In addition, the bottom Yukawa coupling 
increases, and as a consequence m 2 ^ decreases, implying that m A , given by Eq. (3.12), 
also decreases. Indeed, annihilation channels through A exchange are more important 
now and their contributions to the flux will increase it. This effect is shown in Fig. 5, 
which can be compared with Fig. 3. 

The above results concerning the gamma-ray flux are summarized in Fig. 6. There 
we see that the flux for a threshold of 1 GeV is larger than for tan/3 =5. In any case, 
only a very small region of the parameter space turns out to be accessible for GLAST, 
and no region for HESS. It is worth noticing that now only the lower areas bounded by 
the dashed line correspond to a^ USY < 7.1 x 10~ 10 , and therefore would be excluded 
by (g^ — 2) if we impose consistency with the e + e _ data. 

We have also analyzed very large values of tan/3, such as 50. In this case, still 
EGRET data cannot be reproduced, but larger regions of the parameter space would 
be accessible for GLAST and also HESS. In particular, fluxes with 3> 7 < 10 -9 cm -2 
($ 7 < cm~ 2 s" 1 ) can be obtained for a threshold of 1 GeV (60 GeV). 
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Figure 3: Gamma-ray flux <I> 7 for a threshold of 1 GeV, in the parameter space of the mSUGRA scenario 
(to, M) for tan/3 — 5, A = and /i > 0. The dot-dashed curves are contours of <E> 7 in cm" 2 s" 1 units. 
A NFW profile is used with A£l = 10 -3 sr. The region to the left of the light grey (green) dashed line is 
excluded by the lower bound on the Higgs mass. The region to the left of the dotted line is excluded by 
the lower bound on the chargino mass to~± > 103 GeV. The region to the right of the black dashed line 
corresponds to a^ USY < 7.1 x 10~ 10 , and would be excluded by e + e~ data. The region at the bottom (Stau 
LSP) is excluded because the lightest stau is the LSP. The dark grey (red) region fulfils 0.1 < f^o/i 2 < 0.3 
(the black region on top of this indicates the WMAP range 0.094 < f^o/i 2 < 0.129). 



Let us finally remark that the result obtained above concerning EGRET, should 
be expected given the following argument. The relic density can be approximated in 
the simplest situation, when the annihilation cross section is velocity independent, as 
f2/i 2 ~ 3xl0 ~ 2 J ); c ™ 3s ~ 1 . This implies that in order to obtain the result VLh 2 ~ 0.1 one 
needs a cross section av ~ 3 x 10~ 26 cm 3 s . Using this value for the annihilation cross 
section in the galaxy in Eq. (2.3), for a NFW profile with J (l0~ 3 ) w 3 x 10 3 and a 
photon yield of about 10 as can be deduced from Fig. 1, one obtains a gamma-ray 
flux $ 7 ~ 10~ 9 cm~ 2 s _1 , i.e. below EGRET sensitivity. In fact, the value of the 
flux is typically two or three orders of magnitude smaller than this estimation, since in 
the early Universe coannihilations can also contribute to av as discussed above, thus 
3 x 10 _26 cm 3 s _1 can be considered as an upper limit. The results shown in Figs. 4 and 
6 confirm this rough discussion. 
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Figure 4: Scatter plot of the gamma-ray flux $ 7 for a threshold of 1 GeV (left), and 60 GeV (right), as 
a function of the neutralino mass m x , using the same parameter space as in Fig. 3, for tan/3 = 5, A = 
and fi > 0. All points have a^ USY < 7.1 x 10~ 10 , and therefore would be excluded by (g^ — 2) if we impose 
consistency with the e + e~ data. They fulfil all other experimental constraints, and those depicted with light 
grey (magenta) triangles have 0.129 < il^oh 2 < 0.3, with black stars have 0.094 < fi^o/i 2 < 0.129, and 
finally with dark grey (blue) boxes have 0.03 < fi^o/i 2 < 0.094 with the appropriate rescaling of the density 
of neutralinos in the galaxy as discussed in the text. From top to bottom, the first solid line in the left frame 
corresponds to the signal reported by EGRET, and the upper area bounded by the second solid line will be 
analyzed by GLAST experiment (HESS experiment in the right frame). 
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Figure 5: The same as in Fig. 3 but for tan/3=35. The region to the left of the solid line is excluded by 
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4 Departures from mSUGRA 



As discussed in detail in Refs. [54, 55, 56] in the context of direct detection, the 
neutralino-proton scattering cross section can be increased in different ways when the 
structure of mSUGRA for the soft terms is abandoned. Since the diagrams for neu- 
tralino annihilation through neutral Higgs exchange are related to those by crossing 
symmetry, we can use the same arguments here. In particular, it is possible to enhance 
the annihilation channels involving exchange of the CP-odd Higgs, A, by reducing the 
Higgs mass, and also by increasing the Higgsino components of the lightest neutralino. 
As a consequence, the gamma-ray flux will be increased. A brief analysis based on the 
Higgs mass parameters, rn 2 Hd and rnjj^, at the electroweak scale can clearly show how 
these effects can be achieved. 

First, a decrease in the value of the mass of A can be obtained by increasing m 2 H ^ 
(i.e., making it less negative) and/or decreasing rn 2 Hd . This is easily understood from 
Eq. (3.12). 

Second, through the increase in the value of m 2 H an increase in the Higgsino com- 
ponents of the lightest neutralino can also be achieved. Making m 2 H less negative, its 
positive contribution to fj? in (3.10) would be smaller. Eventually \fi\ will be of the 
order of M±, M 2 and x? will then be a mixed Higgsino-gaugino state. Thus annihilation 
channels through Higgs exchange become more important than in mSUGRA, where 

is large and Xi 1S mainly bino. This is also the case for Z-, xf~, and ^-exchange 
channels. 

As pointed out in Refs. [54, 55, 56], non-universal soft parameters can produce the 
above mentioned effects. Let us first consider non-universalities in the scalar masses. 



4.1 Non- universal scalars 

We can parameterise the non-universalities in the Higgs sector, at the GUT scale, as 
follows: 

m 2 Hd = m 2 (l + 5 1 ) , m 2 Hu =m 2 (l + 5 2 ) . (4.13) 

Concerning squarks and sleptons we will assume that the three generations have the 
same mass structure: 

m QL = m 2 (l + 5 3 ), m 2 UR = m 2 (l + <5 4 ) , 



m 



2 m 2 (l + 5 5 ), m 2 dR = m 2 (l + <5 6 ) , 



en 

m\ L = m 2 (l + 5 7 ). (4.14) 



Such a structure avoids potential problems with flavour changing neutral currents 7 . 
Note also that whereas 8{ > —1, i = 3,..., 7, in order to avoid an unbounded from 
below (UFB) direction breaking charge and colour, <5i 5 2 < — 1 is possible as long as 
m 1 = m2 H d + ^ 2 > and ml = m H„ + /W 2 > are fulfilled. 



7 Another possibility would be to assume that the first two generations have the common scalar mass to, 
and that non-universalities are allowed only for the third generation. This would not modify our analysis 
since, as we will see below, only the third generation is relevant in our discussion. 
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An increase in m 2 H at the electroweak scale can be obviously achieved by increasing 
its value at the GUT scale, i.e., with the choice 82 > 0. In addition, this is also produced 
when m,Q L and m 2 lR at Mqut decrease, i.e. taking £3,4 < 0, due to their (negative) 
contribution proportional to the top Yukawa coupling in the renormalization group 
equation (RGE) of rn 2 Hu . Similarly, a decrease in the value of rn 2 Hd at the electroweak 
scale can be obtained by decreasing it at the GUT scale with 8\ < 0. Also, this effect is 
produced when TriQ L and m^ fi at Mqtjt increase, due to their (negative) contribution 
proportional to the bottom Yukawa coupling in the RGE of m 2 H . Thus one can deduce 
that m 2 A will be reduced by choosing also 83^ > 0. In fact non-universality in the Higgs 
sector gives the most important effect, and including the one in the sfermion sector 
the annihilation cross section only increases slightly. Thus in what follows we will take 
Si = 0, i = 3, 7. Let us analyze then three representative cases: 



Clearly, the above discussion about decreasing /j, 2 applies well to case a), where 
the variation in m 2 H through 82 is relevant. This is shown in Fig. 7 for tan (3 = 5 and 
A = 0, which can be compared with Fig. 3. Note that now there is an important area in 
the upper left (no EWSB) where fi 2 becomes negative due to the increasing in 82 with 
respect to the universal case. We see that the decrease of the \x parameter opens a new 
allowed (narrow) region of the parameter space, close to this one. In particular, the 
value of the relic density fi^o is affected due to the increase of the Higgsino components 
of Xi with respect to the dominant bino component of the universal case. Thus the 
relic abundance is placed inside the astrophysical bounds. In Fig. 2 we can see that 
this nature of the LSP enhances the neutralino coupling to the Z (giving rise to bottom 
and top final states when kinematically allowed), and to the lightest chargino \^ an d 
neutralino Xi (inducing W^W^ and ZZ boson final states, especially for low values 



We can see the enhancement in the gamma-ray flux in Fig. 8, where many points 
accessible for GLAST and HESS are obtained. This can be compared with the result of 
Fig. 4 for the universal scenario, where such points are not possible with tan (3 = 5. On 
the other hand, the lower region, below GLAST (and HESS), corresponds to the narrow 
coannihilation branch. In this case the neutralino is mainly bino and the discussion of 
this region for mSUGRA in the previous section applies also here. 

Let us remark, however, that if we impose consistency with the e + e~ data concern- 
ing (g^ — 2), the whole parameter space would be excluded (notice that the region to 
the right of the black dashed line in Fig. 7 corresponds to a^ USY < 7.1 x 10~ 10 ), as in 
the case of mSUGRA. 

For tan (3 = 35, one can see in Figs. 9 and 10 the enhancement of the flux related 
to the coannihilation branch. As discussed in the previous section for mSUGRA, the 
increase is about two orders of magnitude with respect to the case tan (3 = 5 (see Fig. 8), 
as expected by the couplings proportionals to tan/3 in the A-exchange channels. On 
the other hand, one can also see in Fig. 10 that the flux in the Higgsino branch does 



a) 61 

b) 8 1 

c) <5i 



, 82 = 1 , 
-1 , 82 = . 
-1 ,fc = l. 



(4.15) 



of M). 
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not increase significantly. This is because it is mainly produced by the Z exchange that 
does not depend on tan/3 at all. Concerning this branch, it is worth noticing that it is 
not so close to the no EWSB area as in the case of mSUGRA. In particular, one can see 
in Fig. 9 that larger values of the parameters imply a larger gap between the Higgsino 
branch and the no EWSB area. Compare for example the points with m=1000 GeV 
and 2000 GeV, and with the relic density inside the astrophysical bounds. In order 
to obtain the same value of the relic density, for m=2000 GeV a higher value of \i is 
necessary. 

It is worth noticing that now, unlike the tan (3 = 5 case, not all points would be 
excluded by — 2) if we impose consistency with the e + e~ data. In Fig. 10 those to 
the left of the dashed line correspond to a^ USY > 7.1 x 10 -10 , and therefore would be 
allowed on this ground. 

On the other hand, for case b) with tan/3 = 5, studied in Figs. 11 and 12, the only 
region which is allowed by the experimental constraints is the one corresponding to the 

narrow coannihilation branch. Recall that now mh, is not varied and therefore the 

u 

effect of decreasing fj, is not present. As in the case of mSUGRA, the gamma-ray flux 
is small and below the accessible experimental areas. However, for tan [3 = 35, due to 
the effect of decreasing m 2 Hd significantly, m\ can be very small. Notice the area in 
the upper left of Fig. 13 where m 2 A becomes negative. As a consequence, Fig. 14 shows 
similar fluxes to those of the case a) in Fig. 10, but with a larger contribution from the 
A exchange. Thus interesting regions of the parameter space are accessible for future 
experiments. 

Now we can see in Fig. 13 (and Fig. 14) that there are two cosmologically allowed 
corridors. This can be easily understood from the variation in m\. For example, in 
Fig. 13, for a given value of m, and from large to small M, vtia decreases. Starting 
with M = 2000 GeV, one has via > 2m^o, and therefore a small annihilation cross 
section producing a large relic density. However, when tua becomes close to 2m^o the 
annihilation cross section increases and therefore the relic density decreases entering in 
the first allowed corridor with 0.094 < VL^ti 2 < 0.3. For slightly smaller values of M 

Al 

one finally arrives to the A-pole region where tiia ~ 2m^o and the relic density is too 

Al 

small. Smaller values of rriA close this region, producing a decrease in the annihilation 
cross section, which allows to enter in the second corridor with the relic density inside 
the observational bounds. For rriA too small this second corridor is closed, and the 
relic density is too large, ^^oh 2 > 0.3. Finally, one arrives to the region where m 2 A 
becomes negative. 

Concerning case c), this is basically a combination of the previous two cases, a) 
and b). In Fig. 15 one can see that the region where the neutralino is mainly Higgsino 
is already open at tan/3 = 5. This is of course due to the non-universality of rnn u - 
Consequently, large fluxes can be obtained, as shown in Fig. 16. These figures are 
qualitatively similar to those corresponding to case a) with tan/3 = 5, Figs. 7 and 8. 
Moreover, one can see in Fig. 17 and 18 for tan/3 = 35 that the A contribution is now 
also largely present, although the effect of the Higgsino component of the neutralino is 
also important in the region m,M ~ 500 GeV. Notice that the shape of the allowed 
corridors is qualitatively similar to that corresponding to case b), Figs. 13 and 14, 

Let us also mention that the combination of both non-universalities gives rise to 
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very light Higgses, particularly charged Higgses, and that new final states (H + W~ or 
ZH) can be kinematically allowed and even reach 30% of the branching ratio in some 
regions of the parameter space. This effect is not possible in mSUGRA. 

Finally, we have also analyzed very large values of tan/3, such as 50, for cases a), b), 
and c). The results are qualitatively similar, in the sense that many points accessible 
for GLAST and HESS are obtained. 
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Figure 7: The same as in Fig. 3 but for the non-universal scalar case a) Si — 0, S2 — 1, discussed in 
Eq. (4.15). 
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Fi gure 8: The same as in Fig. 4 but for the non- universal scalar case a) Si = 0, #2 = 1, discussed in 
Eq. (4.15). 
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Figure 9: The same as in Fig. 3 but for tan/3 =35 and the non-universal scalar case a) 5\ = 0, 62 = 1, 
discussed in Eq. (4.15). The region to the left of the solid line is excluded by b — * 57. 
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Figure 10: The same as in Fig. 4 but for tan/3 =35 and the non-universal scalar case a) Si = 0, 62 = 1, 
discussed in Eq. (4.15). Now, the region to the right of the dashed line corresponds to a^ u < 7.1 x 10~ 10 , 
and would be excluded by e + e~ data. 
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Figure 11: The same as in Fig. 3 but for the non-universal scalar case b) Si = — 1, 62 = 0, discussed in 
Eq. (4.15). The region to the left of the solid line is excluded by b — > sj. 
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Figure 12: The same as in Fig. 4 but for the non-universal scalar case b) 8\ = — 1, ^ = 0, discussed in 
Eq. (4.15). 



19 





Figure 14: The same as in Fig. 4 but for tan /? = 35 and the non-universal scalar case b) 5\ — —1, 82 = 0, 
discussed in Eq. (4.15). Now, the regions to the right of the dashed lines correspond to a^ USY < 7.1 x 10~ 10 , 
and would be excluded by e + e~ data. These do not include the points with black stars between both dashed 
lines, which are allowed. 
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Figure 15: The same as in Fig. 3 but for the non-universal scalar case c) Si = — 1, 62 = 1, discussed in 
Eq. (4.15). The region to the left of the solid line is excluded by b — > sj. 
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Figure 16: The same as in Fig. 4 but for tan/3 = 35 and the non-universal scalar case c) 81 = — 1, #2 = 1, 
discussed in Eq. (4.15). 
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Figure 18: The same as in Fig. 4 but for tan/3 = 35 and the non-universal scalar case c) 8i = — 1, 82 = 1, 
discussed in Eq. (4.15). Now, the regions to the right of the dashed lines correspond to a^ USY < 7.1 x 10~ 10 , 
and would be excluded by e + e~ data. These do not include the points with black starts between both dashed 
lines, which are allowed. 
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4.2 Non- universal gauginos 

Let us now study the effect of the non- universality in the gaugino masses. We can 
parameterise this as follows: 

Mi = M , M 2 = M(l + 6' 2 ) , M 3 = M(l + <5(,) , (4.16) 

where Mi^.3 are the bino, wino and gluino masses, respectively, and 5^ = corresponds 
to the universal case. 

In order to increase the annihilation cross section it is worth noticing that M3 
appears in the RGEs of squark masses. Thus the contribution of squark masses pro- 
portional to the top Yukawa coupling in the RGE of m 2 H will do this less negative 
if M3 is small. As discussed above, this produces a decrease in the value of | /x | and 
therefore a mixed Higgsino-gaugino lightest neutralino, which may help in order to 
increase the gamma-ray flux. 

Because the mass of the lightest Higgs is very dependent on the value of M3 , through 
radiative corrections, its decrease is very limited. In fact, in order to satisfy the lower 
limit of M3, M in (4.16) may have to increase, thus rather than a decrease in M3 what 
one obtains is an effective increase in Mi and M 2 , which leads to a larger (less negative) 
value of mjj and thus a reduction in the value of |//|. This in turn implies heavier 
neutralinos, when the lightest neutralino is mostly gaugino, and an increase of the 
Higgsino composition, which would be dominant if Mi > \fi\ at the electroweak scale. 
Finally, decreasing the ratio M3/M1 leads to a more efficient neutralino annihilation, 
due to the enhancement in the Higgsino components 8 of x?) entailing a reduction of 

n xS- 

In the following we analyze the case 5' 2 = 0, 6' 3 = —0.4. Lower values of M3 at the 
GUT scale cannot produce a sufficient amount of relic density. We show the situation 
for tan (3 = 5 in Fig. 19. Notice that the Higgs mass bound allows only the region 
corresponding to the narrow coannihilation branch. Thus the A-exchange channel is 
the relevant one here. The fact that in this case the gamma-ray flux is larger than 
in mSUGRA (one can compare Figs. 20 and 4) is due to the different composition 
of the neutralino. For example, for a point with m = 200 GeV and M = 1100 GeV 
(corresponding to Q x =0.16) one obtains Z13 = 0.11 in this non-universal case, whereas 
Z13 = 0.04 in mSUGRA. Thus the Higgsino component of the neutralino is more 
important with non-universality. 

In the large tan (3 regime there is now a new feature. We can see in Fig. 21 that 
we have three cosmologically allowed corridors for tan (3 = 35. To analyze them, let 
us consider a fixed value of m, for example m = 1000 GeV, and go from large to 
small M. For M = 2000 GeV, txia is already smaller than 2m^o (this effect is not 
present for tan/3 = 5), producing a too large relic density. However, for smaller values 
of M the decreasing in ijia is not as steep as in non-universal scalar cases, and for 
M ~ 1500 GeV, tua becomes very close to 2m^o. Thus the relic density is placed 
inside the bounds, 0.094 < ^ x o/i 2 < 0.3. For slightly smaller values of M one arrives to 
the A-pole region since ijia ~ 2m^o , and the relic density is too small. If we continue 

8 Note that the value of m 2 Hd also increases, thus m,A calculated from (3.12) is typically not very affected. 
Although we will see below that the decrease of vtia in some regions can be important for the analysis. 
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to the left, niA > 2m^o, and this close the A-pole region producing a decrease in the 
annihilation cross section, which allows us to enter in the second corridor. Finally, 
the third corridor to the left of Fig. 21 is the usual Higgsino branch generated by the 
increase of the Higgsino components of the neutralino, near the noEWSB area where 
H is small. As discussed in the case of Fig. 9, the gap between this branch and the no 
EWSB area increases with larger values of the parameters. See e.g. m=1500 GeV and 
2000 GeV. It is also worth noticing that now, for a fixed value of m, there are more 
points with a relic density inside the astrophysical bounds. This is so because the value 
of fi is more stable with respect to variations in M. Note in this respect that here we 
are lowering /x through the value of M3, and this is a loop effect. 

Let us finally mention that in the latter region the Z-exchange channel dominates 
producing fluxes with bb final states (tt being kinematically disfavoured). This corre- 
sponds to M smaller than about 900 GeV. The corridors related to the A-pole open 
for M about 600 and 1000 GeV. These three corridors can be clearly distinguished 
in Fig. 22. The mixed Higgsino-gaugino one for 200 < m x < 400 GeV, and the other 
two for 300 < m x < 900 GeV and 450 < m x < 900 GeV. The latter gives rise to the 
largest flux. As in the case of non-universal scalars, we can observe in the figure that 
important regions of the parameter space will be accessible for future experiments. 
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Figure 19: The same as in Fig. 3 but for the non-universal gaugino case 5' 3 
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-0.4, discussed in Eq. (4.16). 
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Figure 20: The same as in Fig. 4 but for the non-universal gaugino case 5' 3 = —0.4, discussed in Eq. (4.16). 
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Figure 21: The same as in Fig. 3 but for tan/3 = 35 and the non-universal gaugino case 5' 3 — —0.4, 
discussed in Eq. (4.16). The region to the left of the solid line is excluded by b — > 37. 
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Figure 22: The same as in Fig. 4 but for tan/3 = 35 and the non-universal gaugino case 6' 3 — —0.4, discussed 
in Eq. (4.16). Now, the regions to the right (and below) the dashed lines correspond to a^ USY < 7.1 x 10~ 10 , 
and would be excluded by e + e~ data. These do not include the points with black stars between both dashed 
lines, which are allowed. 
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4.3 General case: non-universal scalars and gauginos 

Let us now consider the general case where the soft super symmetry-breaking terms for 
both scalar and gauginos have a non-universal structure. This was studied in detail 
for the case of direct detection in Ref. [56], and we will follow closely the analysis. 
We will be mostly interested in analyzing possible departures with respect to the pa- 
rameter space which was allowed in the previos subsections, producing large values of 
the neutralino annihilation cross section. For this reason, we will concentrate on the 
interesting choices for scalar non-universalities given in Eq. (4.15), and study the effect 
of adding gaugino non-universalities to these. 

We will consider the possibility of increasing the value of M3 with respect to Mi, 
which can be done with 5' 3 > in Eq. (4.16). In this case, the constraint on the Higgs 
mass and on b — > 57 will be satisfied for smaller values of M, and therefore the effective 
value of Mi can be smaller than in the universal case. Thus lighter neutralinos can be 
obtained. 

In particular, we will study the theoretical predictions for the gamma-ray flux for 
an example with 5' 2 3 = 3, tan f3 = 35, A = 0, and the three choices of Higgs non- 
universalities of Eq. (4.15). As we will see, this choice of gaugino parameters favours 
the appearance of light neutralinos within the GLAST sensitivity range. In these cases 
the typical values for the relic density are too large and therefore not consistent with 
the WMAP result. A very effective decrease can be achieved when the h and Z-poles 
are crossed at m^o = rrih/2 and Mz/2, respectively, giving rise to a very effective 
neutralino annihilation through the corresponding s-channels. This is evidenced by 
the narrow chimneys in the cosmologically preferred regions in Fig. 23, where case c) 
of Eq. (4.15) is shown. Because of the very efficient decrease in the CP-odd Higgs mass 
in this case, annihilation of very light neutralinos can be boosted and thus the correct 
relic density obtained. This allows the existence of neutralinos with m^o ~ 30 GeV 
which produce a flux even close to EGRET sensitivity, as shown in Fig. 24. For the 
other two cases, a) and b) of Eq. (4.15), the figures are qualitatively similar, although 
the fluxes are slightly smaller, bounded by 10~ 9 cm~ 2 s _1 . For the three cases these 
light neutralinos can in principle be detected by GLAST. Of course, they cannot be 
detected by atmospheric Cherenkov telescopes, as e.g. HESS, since these have a large 
threshold energy of about 50 GeV. 

It is worth noticing that the existence of very light staus with light neutralinos in 
case a) with ^,3 = 3, makes it possible the domination of the lightest stau exchange 
channel with respect to the A one (see Fig. 2) in regions of the parameter space, 
producing a large amount of taus. This comes from the RGE governing the mj, 3 
evolution. Indeed, decreasing Mi for a given M2 decreases the value of rriL 3 at the 
electroweak scale, giving a lighter stau. On the other hand, cases b) and c) do not 
exhibit such a behavior because decreasing at the same time m# d cancels the effect 
(the two contributions appear with an opposite sign in the RGE). 

Let us now discuss the possibility of obtaining even lighter neutralinos. One of the 
requirements for the appearance of such very low neutralinos is to have Mi <C /j,, M2 at 
low energy (thus having almost pure binos). This can be achieved with adequate choices 
oibb gaugino non- universalities, in particular with 5' 23 3> 1. However, as mentioned 
above, without a very effective reduction of tua, the relic density would be too large, 
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and therefore inconsistent with observations. Here the presence of non- universal scalars 
is crucial. In particular, non-universalities as the ones we have described in the Higgs 
sector in Eq. (4.15) provide a very effective way of lowering tjia and are thus optimal 
for this purpose. 

More specifically, it is in case b) and especially c) where the reduction in iua is more 
effective (not being so constrained by regions with fi 2 < 0) and for this reason very 
light neutralinos can easily appear. We have already seen in the previous example how 
this happens for case c) with tan (3 = 35 and 5' 2 3 = 3 (see Fig. (24)). On the contrary, 
in case a) higher values of tan (3 are required in order to further reduce the value of m^. 
One can check explicitly that tan/3 > 33 is sufficient to obtain m^o < Mz/2 in cases 
b) and c), whereas tan/3 > 45 is necessary in case a). In all the three cases 5' 23 > 3 
leads to these results. 

Let us therefore complete our discussion by analysing the case 5' 2 3 = 10 for the 
choice a) of scalar non-universalities in Eq. (4.15). The result is shown in Fig. 26 for 
tan (3 = 50. We observe the appearance of very light neutralinos, whose cross section 
can be in range of detectability of near-future experiments. In particular, points with 
a gamma-ray flux in the region of GLAST sensitivity are obtained with m^o ~ 15 
GeV. As mentioned for 5' 2 3 = 3, these points cannot be detected by HESS because 
of the large threshold energy of about 50 GeV. For the other two cases, b) and c) of 
Eq. (4.15), the figures are qualitatively similar. 

The effect of the different constraints on the corresponding (m, M) parameter space 
is represented in Fig. 25. Note that due to the large value of 6' 3 the resonances with 
the lightest Higgs and the Z are not sufficient to place the relic abundance inside the 
bounds, and the regions giving rise to very light neutralinos with a consistent relic 
density are extremely narrow. In these points the mass of the CP-odd Higgs can 
be very close to its experimental limit, uia ^ 100 GeV. Since the neutralino mass is 
so small, the region excluded due to the neutralino not being the LSP is negligible. 
However, now the lower bound on the stau mass plays an important role. In fact, an 
important region of the parameter space is excluded for having m~ < 0. 
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M (GeV) 

Figure 23: The same as in Fig. 3 but for tan/3 =35 and the non-universal scalar and gaugino case 
82 = 1) Si = —1, 5' 2 = 8' 3 = 3. The region to the left of the solid line is excluded by b — > S7. 
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Figure 24: The same as in Fig. 4 but for tan/3 =35 and the non-universal scalar and gaugino case 
82 = 1,8\ = — 1, <5 2 — ^3 — 3, using the same parameter space as in Fig. 23. Now, the region to the right of 
the dashed line corresponds to a^ USY < 7.1 x 10~ 10 , and would be excluded by e + e~ data. 
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Figure 26: The same as in Fig. 4 but for tan/3 = 50 and the non-universal scalar and gaugino case 



Mi = 0, 8' 2 



6', 



10, using the same parameter space as in Fig. 25. Now, in the left frame, the 



region to the right of the dashed line corresponds to a^ USY < 7.1 x 10 10 , and would be excluded by e + e 
data. In the right frame all points would be excluded by this constraint. 
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5 Conclusions 



We have analyzed the theoretical predictions for the indirect detection of neutralino 
dark matter in the context of a general SUGRA theory, where both the scalar and 
gaugino soft supersymmetry-breaking terms can have a non-universal structure. In 
particular, we have computed the predictions for the gamma-ray fluxes arising from 
the galactic center due to neutralino annihilation, and compared it with the sensitivity 
of detectors. 

Although the presence of non-universal scalar and gaugino masses is not able to 
explain the intriguing signal detected by the space-based experiment EGRET, the an- 
nihilation cross section increases significantly with respect to the universal (mSUGRA) 
case, producing larger gamma-ray fluxes. For example, whereas mSUGRA is not com- 
patible for tan (3 < 35 with the sensitivity of another crucial projected experiment, 
GLAST, non-universality in the scalar sector can increase the gamma-ray fluxes about 
two orders of magnitude, entering in the region which will be analyzed by this ex- 
periment. For the atmospheric Cherenkov detector HESS, which has already begun 
operations, the situation is qualitatively similar. Of course, for larger values of tan (5 
a similar increase in the gamma-ray fluxes is obtained. The above effects are mainly 
due to the Higgsino nature of the neutralino and / or the decrease of the CP-odd Higgs 
(A) mass in some regions of the parameter space with non- universalities. 

When non-universal gauginos are considered, one can also obtain regions of the 
parameter space with annihilation cross sections much larger than in mSUGRA. This 
is also the case when it is combined with non-universal scalars, i.e. the most general case 
in SUGRA. In addition, in this general case, departures from the allowed parameter 
space obtained in the previous situations are possible. For example, neutralinos as 
light as m^o ~ 15 GeV, and with large fluxes can be obtained. 

This general analysis can be very useful in the study of more specific cases, such 
as the SUGRA theories resulting at the low energy limit of string constructions. In 
particular, D-brane scenarios in type I string theory give rise to theories where non- 
universalities appear both in the scalar and gaugino sectors. 

Let us finally remark that the above results have been obtained using the standard 
NFW density profile for our galaxy. However, we have discussed how they can easily 
be extended for other widely used profiles, such as Isothermal, Kravtsov et al., Navarro 
et al., and Moore et al. In particular, one can use Table 1 where the function J(A$7) 
has been computed for the different profiles (and possible angular resolutions). This 
simply enters as a scaling factor in the gamma-ray fluxes. Obviously, a Moore et al. 
profile, which is about two orders of magnitude larger than the NFW one, would modify 
drastically the results concerning EGRET. 
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